Reply to Comment on "Material evidence of a 38 MeV boson' 
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We reply to a very recent Comment [3| by Bernhard, Friedrich, Schliiter, and Schonning for the 
COMPASS Collaboration, in which it is stated that Monte-Carlo simulations of COMPASS data do 
not support any interpretation of a peaked structure in diphoton invariant-mass distributions below 
the Tv'^ mass in terms of a new resonance. Here we show, by directly comparing the simulations to 
the COMPASS data themselves, that the authors' claim is not substantiated by the Monte-Carlo 
results presented in the Comment. 

PACS numbers: 11. 15. Ex, 12.38.Aw, 12.39.Mk, 14.80.Ec 



Before directly replying to the Comment [ij of Bern- 
hard, Friedrich, Schliiter, and Schonning (BFSS) (for the 
COMPASS Collaboration) on our interpretation [2|pf a 
resonance-like structure in recent COMPASS data [3, 3 1 
let us first recapitulate some of the essential features of 
our work, for those readers who are not so familiar with 
it. 



MOTIVATION 

In Ref. an iS'0(4, 2) conformally symmetric model 
was proposed for strong interactions at low energies, 
based on the observation that confinement can be 
described by an anti-De Sitter (aDS) background geom- 
etry. The possibility of such a strategy had already been 
studied, almost a century ago, by H. Weyl 0], who found 
that the dynamical equations of gauge theories retain 
their flat-space-time form when subject to a conformally- 
flat metrical field instead of the usual Minkowski back- 
ground. The unification of electromagnetism and strong 
interactions can be justified by the very subtle balance 
between these forces in the nucleus, where just one neu- 
tron more or less can make the difference between stabil- 
ity or instability. 

Confinement of quarks and gluons is modeled by the 
introduction of two scalar fields which spontaneously 
break the SO{A, 2) symmetry down to SO{3, 2) and 
50(3)(X)S'0(2) symmetry, respectively. Moreover, asym- 
metric second-order tensor field is defined that serves as 
the metric for flat space-time, coupling to electromag- 
netism. Quarks and gluons, which to lowest order do not 
couple to this tensor field, are confined to an aDS uni- 
verse [l^, having a finite radius in the flat space-time. 
This way, the model describes quarks and gluons, which 
oscillate with a universal frequency — independent of the 
flavor mass — inside a closed universe, as well as photons, 
which freely travel through flat space-time. 

The flelds in the model of Ref. Q comprise one real 



scalar fleld a and one complex scalar fleld A. Their dy- 
namical equations were solved in Ref. Q for the case 
that the respective vacuum expectation values, given by 
(jQ and Aq, satisfy the relation 



kol > |Ao| 



(1) 



A solution for ctq of particular interest leads to aDS con- 
finement, via the associated conformally flat metric given 
by cnjt^i'- 

The only quadratic term in the Lagrangian of Ref. [5| 
is proportional to 



(2) 



Hence, under the condition of relation ([T]), one obtains, 
after choosing vacuum expectation values, a light a field 
associated with confinement, and a very heavy complex 
A field associated with electromagnetism. Here, we will 
study the — supposedly light — mass of the scalar field 
that gives rise to confinement. 

The conformally symmetric model of Ref. Q in itself 
does not easily allow for interactions between hadrons, 
as each hadron is described by a closed universe. There- 
fore, in order to compare the properties of this model to 
the actually measured cross sections and branching ra- 
tios, the model has been further simplified, such that only 
its main property survives, namely its flavor-independent 
oscillations. This way the full aDS spectrum is, via light- 
quark-pair creation, coupled to the channels of two — 
or more — hadronic decay products for which scatter- 
ing amplitudes can be measured, thus relating observed 
resonances to the spectrum of aDS. 

The aDS spectrum reveals itself through the structures 
observed in hadronic invariant-mass distributions. How- 



ever, as we have shown in the past (sec Ref. and ref- 
erences therein), there exists no simple relation between 
enhancements in the experimental cross sections and the 
aDS spectrum. Nevertheless, this was studied in parallel, 
for mesons, in a coupled-channel model in which quarks 
are conflned by a flavor-independent harmonic oscillator 
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12lll3l|. Empirically, based on numerous data on mesonic 
resonances measured by a large variety of experimental 
collaborations, it was found 1J| that an aDS oscillation 
frequency of 



190 MeV 



(3) 



agrees well with the observed results for meson-meson 
scattering and meson-pair production in the light [l5| . 
heavy- light [31, and heavy 17, [3| flavor sectors, thus 
reinforcing the strategy proposed in Ref. Q. 

A further ingredient of the model for the description 
of non-exotic quarkonia, namely the coupling of quark- 
antiquark coniponents to real and virtual two-meson de- 
cay channels [l9| via ^Pq quark-pair creation, gives us a 
clue about the size of the mass of the a field. For such 
a coupling it was found that the average radius tq for 
light-quark-pair creation in quarkonia can be described 
by a flavor-independent mass scale, given by 



(4) 



where fi is the effective reduced quarkonium mass. In 
earlier work, the value po ~ y^JEovQ = 0.56 12, 1^ was 
used, which results in M = 30 MeV for the correspond- 
ing mass scale. However, the quarkonium spectrum is 
not very sensitive to the precise value of the radius tq , in 
contrast with the resonance widths. In more recent work 
[20I 21 1, slightly larger transition radii have been applied, 
corresponding to values around 40 MeV for M. Never- 
theless, values of 30-40 MeV for the flavor-independent 
mass M do not seem to bear any relation to an observed 
quantity for strong interactions. However, in Refs. [^.[22^ 
we have presented experimental evidence for the possible 
existence of a quantum with a mass of about 38 MeV, 
which in light of its relation to the '^Pq mechanism we sup- 
pose to mediate quark-pair creation. Moreover, its scalar 
properties make it a perfect candidate for the quantum 
associated with the above-discussed scalar field for con- 
finement. 



distributions of e+e" 
signal 
claims. 
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and finally hybrid 



III and 



has much too low statistics to make firm 
However, we noted that all these results point 
in the same direction. 

Indeed, it seems highly unlikely that by sheer coici- 
dence one finds the same kind of oscillations in four dif- 
ferent sets of data 0, [i^l involving different fiavors, sta- 
tistical fluctuations at ±38 MeV in yet another four sets 
of different data, and finally a resonance-like fluctuation 
at 38 MeV in a further set of data. Moreover, the re- 
sulting scalar mass comes out where it was predicted via 
our analyses in mesonic spectroscopy (see Ref. 
references therein). 



ENHANCEMENT AROUND 40 MEV IN 
COMPASS DATA 

In Ref. 0, we presented more pieces of evidence, one 
of which, viz. a small resonance-like signal in the dipho- 
ton invariant-mass distribution Q of the COMPASS Col- 
laboration, is considerably more conclusive than previ- 
ously reported [22| signals, owing to much higher statis- 
tics. The data, depicted in Fig. [1] have been obtained 
at the two-stage magnetic COMPASS spectrometer at- 
tached to the SPS accelerator facility at CERN These 
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FIRST OBSERVATIONS 

In recent papers 0, [2^ , we have presented a variety 
of indications of the possible existence of a light boson 
with a mass of about 38 MeV, henceforth referred to 
as E{38). These indications amounted to a series of 
low-statistics observations all pointing in the same di- 
rection. Each of the results alone, viz. interference ef- 



fect [23[ , small flavor- independent oscillations in electron- 
positron and proton-antiproton annihilation data (j^ . 
the excess signals visible [2 
butions of T (2 35i) ^ 7r+^ 
T{3^Si) -> 7r+7r-T(l3S'i 



|22l in the 

-T(l3^l) ^71 



TT+TT-T (2^Si) n+TT' 



mass distri- 
fi^ fi~ , in 
fi^ fi^ , in T (3351) 
and in the e^e~ mass 



FIG. 1: Detail of the invariant two-photon mass distribution 
of COMPASS [3. We observe a clear signal, with maximum 
at about 39 MeV. 

data seem to have enough statistics to substantiate the 
existence of a light boson with mass around 40 MeV. 
However, in a more recent version 0] of their work, the 
COMPASS Collaboration added the following remark to 
the figure caption regarding the enhancement in Fig. [TJ 

"The structures below the tt*^ mass peak are 
artefacts of low energetic photon reconstruc- 
tion due to secondary interactions in the de- 
tector material and to cuts in the reconstruc- 
tion algorithm. They should not be mistaken 
for any physical signal. " 
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Although it may of course be possible to obtain 
resonance-like structures by artefacts in data collection, 
we are convinced this is not the case for the signal in 
the 40 MeV region, because it coincides surprisingly well 
with the other observations reported in Refs. [23|. 

Furthermore, it is clear that the light boson cannot 
be an ordinary hadron emerging from a hadronic vertex, 
unless at an extremely low rate. Otherwise, it would 
have been observed long ago. Ordinary hadronization 
in high-energy collisions gives rise to pions, kaons, and 
other hadrons that are stable with respect to strong in- 
teractions. These are processes in which quark-pair cre- 
ation and gluonic jets dominate. But on the other hand, 
judging from the amount of events in the low-mass en- 
hancement in the two-photon data, which is about 10% 
of the number of events in the rj signal, it does seem to be 
produced with a reasonably large rate in the COMPASS 
experiment. Such a rate indicates that it most probably 
is a hadronic particle, though with very peculiar proper- 
ties that still have to be understood. 

In this paper we shall assume that the £'(38) boson has 
the shape of a spherical bubble, as predicted by anti-Dc 
Sitter confinement H,!!^, i.e., a thin film of glue. Due to 
possible surface oscillations, such a system almost has a 
continuum of excitations. If we assume that the stability 
of the bubble rapidly decreases for excitations and allow 
for a spreading of some 20 MeV, then we may fully re- 
construct the line shape shown in Fig. [TJ Moreover, we 
point out that the small excess in the diphoton data ob- 
tained in a different experiment [2^ also satisfies the here 
proposed distribution, albeit with much lower statistics. 

A further important ingredient for the reconstruction 
of the present event distribution is the experimental res- 
olution. Given the 77 signal in the COMPASS ^ data, 
which is much broader than that of tt", and the fact that 
7r° and 77 have decay widths considerably smaller than 
the bin size of 0.5 MeV, we may assume that the experi- 
mental resolution is better at lower invariant mass. 

In order to proceed, we compare in Fig. [2] the three 
observed relevant structures in the COMPASS diphoton 
data, representing the TfO, 77, and £(38) bosons. The 
three enhancements are scaled in width. Also, the very 
different heights have been adjusted in Fig.[2l so as to al- 
low comparing the line shapes. The data spreading must 
be entirely due to the experimental resolution. Further- 
more, we assume that the £(38) boson is reasonably sta- 
ble and so should have a spreading comparable to the 
other two enhancements. This seems to be confirmed by 
Fig. [21 except for the long £(38) tail at higher invariant 
mass. 

From the representation of the enhancements in Fig.[2l 
we infer an experimental resolution of about 5.2 MeV 
in the 40 MeV region. The convolution of the two 
Gaussians, one for the spreading in the masses of the 
light bosons and another for the experimental resolution 
around 40 MeV, for 56,000 generated Monte Carlo (MC) 
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FIG. 2: The three enhancements observed by us in the 77 
invariant mass distribution of Ref. tt" (red), 77 (green), and 
£(38) (thin black curve). The relative heights and widths of 
the three signals are in reality different, viz. tt" : 77 : E{3S) ~ 
1 : 0.06 : 0.03 and 1 : 3.1 : 0.5, respectively, but have been 
adjusted here in order to compare the line shapes. 



events is show in Fig. [3] 

The here followed procedure is extremely simple, but it 
reproduces the experimental result quite well (see Fig. [3]). 
In particular, it describes how the peak shifts to a lit- 
tle bit higher value than the expected 38 MeV, and fur- 
thermore reproduces the observed tail for larger invariant 
masses. The central mass of the light boson allows for 
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FIG. 3: Invariant two-photon mass distribution of Ref. 0] 
(green), compared to model result (black). 

variations of a few tenths of MeV, from which we may 
estimate its mass to lie between 38.0 and 38.4 MeV. 

We may thus conclude, with very simple assumptions, 
that the event distribution of COMPASS near 40 MeV 
may be understood by the production of a light boson 
with a mass of 38.2±0.2 MeV. 

The DELPHI Collaboration performed an analysis [2^ 
of inclusive tt*^ production in Z'^ decays. The diphoton 
decays of 7r° were reconstructed by using pairs of com- 
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binations of converted photons as well as HPC photons, 
which are photons that were reconstructed in the barrel 
electromagnetic calorimeter of the DELPHI detector. In 
particular, reasonably high-statistics diphoton data were 
shown by the DELPHI Collaboration in Fig. 3b (1 con- 
verted photon and 1 HPC photon) of Ref. [2^, which 
figure also contains a rather precise background fit. 

We apply a similar strategy to the residual signal of the 
DELPHI data as we did above to the COMPASS data of 
Ref. (3]. The resuh of the DELPHI data with an 
assumed resolution of 10 McV, is depicted in Fig. So 
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FIG. 4: Invariant two-photon mass distribution (blue) and 
the background curve (black) of Ref. [2^ . Bottom: difference 
between the data and the DELPHI background curve, multi- 
plied by 5. Thus, 1000 at the vertical axis represents 200 for 
the dots and error bars (red). The model result for 1400 MC 
data, represented by the lower histogram, has been shifted 
downwards by 51 events/5 MeV, in order to compensate for 
a mismatch between the background curve and the data at 
lower invariant mass. 

it seems to us that our assumption on the existence of 
a scalar boson with a mass of about 38.2 MeV, is very 
plausible, since the setup and experimental conditions of 
the DELPHI experiment are very different from those of 
COMPASS. 

SIMULATION OF DIPHOTON DATA 

In principle, the diphoton invariant-mass distribution 
below the nominal tt*^ mass has a structure as represented 
by the red curve in Fig. [SJ However, the two electromag- 
netic calorimeters, ECALl at about 6 meters from the 
target and ECAL2 at about 30 meters downstream, do 
not accept low-energy photons. This results in the non- 
observation of the enhancement at zero invariant dipho- 
ton mass. Hence, at very low masses no events are ob- 
served, since they are removed by the trigger system of 
the EM calorimeters, as indicated by the yellow area in 
Fig. [S] It has the effect that at low mass a peak shows 
up in the data. 



In FigJ5]we depict the low-mass peak for the 7r~p data 
of Ref. [3[, while in the inset of Fig. [SJ we show the low- 
mass peak for the pp data of Ref. [J] . Note that these low- 
mass peaks are very different for different experiments, 
namely the one in the pp data of Ref. Q is roughly 3.8 
times larger than the signal at 25 MeV, whereas the one 
in the Tr~p data of Ref. [3| is only about 1.4 times higher 
than the signal at 25 MeV. 
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FIG. 5: Invariant two-photon mass distributions below the 
nominal tt" mass. The red curve in principle indicates the 
general aspect of such a distribution. Data removed by the 
trigger system are represented by the yellow area. So the 
green area in principle represents the final data. The resulting 
low-mass peak of the main figure coincides with the tt^p data 
of Ref. [3], whereas, in the inset, we display the low- mass peak 
for the pp data of Ref. Qj. 

Now, the data selection system does of course also in- 
fluence the aspect of the data for higher diphoton masses. 
It may even result in some structures which resemble 
resonances. Moreover, several physical processes in the 
experimental setup, not related to the purpose of the ex- 
periment, may result in further structures. 

The authors of Ref. P gave the following list of mech- 
anisms that may result in structures in their data. 

• Secondary tt^ mesons produced in the detector ma- 
terial downstream of the target lead to diphoton 
masses which arc below the nominal tt'^ mass when 
reconstructed assuming a target vertex. 

• Material concentrated in detector groups leads to 
peak-like structures. 

• Secondary e^e^ pairs from photon conversion in 
the spectrometer material lead to low-mass struc- 
tures. 
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• Cuts applied in the reconstruction software lead to 
additional structure in the low-mass range. 

Most of those processes occur, of course, in the EM 
calorimeter ECAL2, because it is further downstream 
than ECALl, and therefore picks up more contamina- 
tion due to unwanted processes. 

These artefacts are reproduced in the MC simulation 
of Ref. [l| for the reactions under study, using a com- 
plete description of the spectrometer material and em- 
ploying the same reconstruction software as for the real 
data analysis. 

In Fig. iniwe reproduce the three independent simula- 
tions reported in Ref. We also add up these three 
simulations, with relative weights ECALl : ECAL2 : 
(ECAL1-HECAL2) = 0.48 : 1.00 : 0.86 for both 7s in 
ECALl, both 7s in ECAL2, and one 7 in each of the 
two electromagnetic calorimeters, respectively, in order 
to have an idea how well they represent the true data 
of Ref. 0. It should be noted, however, that summing 
up three independent simulations is not the same as per- 
forming an MC simulation for the complete experimental 
setup. Therefore, our procedure is not totally reliable, 
though it is the best we can do with the available infor- 
mation. 
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FIG. 6: The three diphoton MC distributions of Ref. [l[ simu- 
lating low-energy 77 production, with different weight factors. 
Black curves: both 7s in calorimeter ECALl or one in ECALl 
and one in ECAL2; green: both 7s in ECAL2. Weight factors: 
ECALl : ECAL2 : (ECAL1+ECAL2) = 0.48 : 1.00 : 0.86. 
Blue curve: sum of green and two black curves. Red curve: 
the experimental data of Ref. 01 . 

We observe from Fig. IH] that for some reason the tt" 
signal, which in the experimental data sets out at about 



50 MeV, is not included in the simulations [ij. Because 
of this omission, the simulated data for one 7 observed 
in ECALl and another 7 in ECAL2 have to contribute 
more than what one would expect from the COMPASS 
setup. 

In Fig. [7] we do include the experimental signal of the 
7r°, with its mass at 135 MeV and an average resolution 
of 12.7 MeV. Once again, this is of course not the cor- 
rect procedure. But it is all we have at our disposal for 
a reconstruction of the experimental data. The simula- 
tions of Ref. [H are here taken in proportions ECALl : 
ECAL2 : (ECAL1-HECAL2) = 0.80 : 0.93 : 0.25 of the 
reported MC data, which seem more reasonable than the 
previous weights employed in Fig. [6l 
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FIG. 7: As FIG. H but now with weight factors ECALl : 
ECAL2 : (ECAL1+ECAL2) = 0.80 : 0.93 : 0.25, and a simu- 
lation of the tt" signal added in the blue curve. 

At this point, we may conclude that experiment is well 
reproduced for diphoton invariant masses larger than 50 
MeV. Nevertheless, we wonder how the MC distributions 
of Ref. [l[ would work out for higher masses, up to 1.0 
GeV. However, for masses below 50 MeV, the simulations 
of Ref. [H do not at all agree with experiment. 

In particular, the left-hand peak, near an invariant 
diphoton mass of 10 MeV, is not reproduced at all. This 
signal, which is almost entirely due to data collected at 
the EM calorimeter ECAL2, depends rather delicately on 
cuts in the data selection, as we explained in the forego- 
ing. The result of Ref. [l[ leaves us with the impression 
that, in order to specifically enhance selected artefacts 
from ECAL2 in the 40 MeV region, the applied cuts in 
the MC simulation are exagerated with respect to those 
actually employed in experiment. However, it is also pos- 
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sible that the artefacts produced in the COMPASS setup 
are just not well studied yet for small diphoton masses. 

The MC simulations shown in Ref. [l| were actually 
produced for the process studied in Ref. However, 
in that study the background-to-signal ratio for the en- 
hancement near 40 MeV is the same as for the process 
considered in Ref. Q. For completeness, we include in 
Fig. El a comparison of the MC simulations of Ref. [ij 
and the actual experimental data of Ref. We find 
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FIG. 8: Similar to FIG. [7 but now compared to the actual 
experimental data of Ref. J . 

that also in this case the experimental data below 50 
MeV arc not at all well described by the MC simulations 
of Ref. m . The discrepancy between data and simulation 
at low diphoton masses is even more serious. 

In conclusion, we welcome the efforts of the authors of 
Ref. [l| to explain the observed enhancement in the 35- 
51 MeV diphoton invariant-mass region by conventional 
methods. However, with the present MC simulations, 
the existence of a resonance-like structure at about 38 
MeV cannot be excluded at all. We suggest to include 
in future simulations the possibility of £'(38) dcsintegra- 
tion in the EM calorimeter ECAL2, since we expect that 
this tentative light boson has rather stable modes, which 
could easily survive the 30 meters that separate ECAL2 
from the target. 



FUTURE 

Although the question whether there exists a (scalar) 
boson with a mass of about 38 MeV does not depend ex- 



clusively on the existence of a resonance-like structure in 
the experimental data of Ref. Q , it is at present the clear- 
est signal we have found in many experiments. Diphoton 
data for the mass interval 10-100 MeV are very rare and 
usually with low statistics. Therefore, it is very impor- 
tant that the present issue be settled, which requires a 
profound understanding of all possible sources of arte- 
facts. 
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